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We introduce a novel approach to scale-decomposition of the fluid kinetic energy 
(or other quadratic integrals) into band-pass contributions from a series of length- 
scales. Our decomposition is based on a multiscale generalization of the "Germano 
identity" for smooth, graded filter kernels. We employ this method to derive a 
budget equation that describes the transfers of turbulent kinetic energy both in 
space and in scale. It is shown that the inter-scale energy transfer is dominated 
by local triadic interactions, assuming only the scaling properties expected in a 
turbulent inertial-range. We derive rigorous upper bounds on the contributions 
of non-local triads, extending the work of Eyink (2005)1^ for low-pass filtering. 
We also propose a physical explanation of the differing exponents for our rigorous 
upper bounds and for the scaling predictions of Kraichnan (1966,1971)'2E1. The 
faster decay predicted by Kraichnan is argued to be the consequence of additional 
cancellations in the signed contributions to transfer from non-local triads, after 
averaging over space. This picture is supported by data from a 512'^ pseudospectral 
simulation of Navier-Stokes turbulence with phase-shift dealiasing. 
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I Introduction 



The traditional Richardson-Kolmogorov-Onsager picture of turbulence supposes an en- 
ergy cascade local in scale. That is, the triad interactions responsible for energy flux 
across a given length-scale are supposed to involve three modes all of comparable scale, 
differing at most by some fixed ratio. This statement also implies that the energy transfer 
is between modes at roughly the same scale. The scale-locality property is fundamental 
to justify the concept of universality of small-scale statistics in turbulent flow. In a re- 
cent wot)^ one of us has given a rigorous proof of the scale-locality property for solutions 
of the Navier-Stokes equation, based on the sole assumption that the velocity field has 
the regularity/scaling properties observed in turbulent flow. The argument employed a 
smooth coarse-graining (or low-pass filtering) technique to resolve the velocity fields si- 
multaneously both in space and in scale. The results in Eyink (2005)'^ apply to individual 
Navier-Stokes solutions, with no statistical averaging, pointwise in space. They therefore 
give substantial support to the classical picture of a local energy cascade. 

Recent numerical studies^ ^ * ^ * ^ * ^ * ^ ! have largely supported these conclusions about en- 
ergy flux. All of these groups verify the individual dominance of nonlocal triadic inter- 
actions but conclude that the net energy flux from summation over triads is dominated 
by local interactions. In particular, the studies'^'^ calculated the locality function W{s) 
of Kraichnanl^ which measures the contribution to energy flux from wavenumber triads 
{k,q,p) with a fixed scale-ratio s = m.m{k,q,p}/med{k,q,p}. They verify to reasonable 
accuracy the prediction of Kraichnan that W{s) oc s^/^ for s <^ 1, in agreement with 
some earlier DNS results'^^'^. Note that this decay of nonlocal contributions in s is even 
faster than that proved in Eyink (2005)'^, which gives W{s) = 0{s'^^^) as a rigourous 
upper bound. The studies'^'^ find also no difference between graded and sharp filters, 
at least with bands defined for a geometric sequences of scales. 

Nevertheless, debate continues on this important issue. In addition to using low-pass 
filtering as in Eyink (2005)'^, the numerical works^ ^ * ^ * ^ * ^ ! have studied the problem also 
by using band-pass filtering. These studies employed both sharp-spectral and graded 
kernels to define modal bands. The relevant nonlinear transfer functions for band-pass 
fields are no longer energy flux but interband energy transfers. In this setting, one 
gj^Q^pHMl j^g^g claimed to verify the importance of nonlocal triadic interactions. Their 
DNS results at resolutions up to 2048^ support earlier observations of "local transfer by 
nonlocal triads" I^^IIIIIIIIS, This really should be called ultralocal transfer because it is 
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between wavenumber modes k,p differing by a fixed amount p — k = ko, where ko ^ 1/L 
is the low wavenumber at the energy peak. This is not to be confused with the traditional 
idea of local transfer, in which the wavenumber ratio p/k is fixed and of order 2. 

How does one explain the apparent discrepancy between these various DNS results 
and the theoretical analysis in Eyink (2005)^? At least three possibilities come to mind. 
One is that there may be an essential difference between transfers for band-pass filtered 
fields and energy flux for low-pass filtered fields. A second possibility is that the key 
difference is between sharp spectral and graded filters. Finally, a third possibility is that 
the important difference is between wavenumber bands of constant size on a linear scale 
versus those with constant size on a logarithmic scale. We stress that the rigorous results 
in Eyink (2005)'^ were proved for graded low-pass filters with logarithmic binning (and 
logarithmic grading of the filter kernel). The logarithmic wavenumber scale is necessary 
in order to achieve simultaneous resolution of the dynamics both in space and in scale. 

In this paper we shall address the first possibility, by extending the rigorous results 
of Eyink (2005)12 to a band-pass decomposition via dilated, smooth filters. As in the 
previous low-pass analysis, we shall show that there are simple equations that describe 
the dynamics of the band-pass energy densities, all terms of which have intuitive physical 
interpretations. Furthermore, we show that band-to-band transfers are local in both scale 
and in space by deriving rigorous upper bounds on the nonlocal triadic contributions. We 
shall also explain on physical grounds the difference between the rigorous upper bounds 
and the more refined scaling predictions of Kraichnanf^. In essence, the contribution of 
the nonlocal triads is smaller than implied by those bounds due to decorrelation of dis- 
parate scales, which implies substantial cancellations of signed contributions to transfer. 
This idea was already proposed by Tennekes & Lumley'^ (Section 8.2) and Kraichnan 
(1971)"^. We shall demonstrate these cancellation effects with data from a 512'^ DNS 
of Navier-Stokes turbulence. These results taken all together confirm the conclusions of 
Eyink (2005) ^ that turbulent energy cascade is a scale-local process, in agreement with 
the traditional view. Furthermore, they show that analysis of band-pass filtered ener- 
gies give results essentially the same as the analysis of low-pass filtered energies. In the 
second part of this study'^ (hereafter referred to as II) we shall consider the remaining 
possibilties, related to sharp-spectral versus graded filters and linear versus logarithmic 
scale-resolution. These turn out to be more subtle and difficult issues to explicate. 

The outline of the paper is as follows. In section II we present our coarse-graining 
approach. Section II. A gathers some material on low-pass filtered fields and their energy 
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budgets. Section II. B applies these techniques to develop a band-pass decomposition 
of the kinetic energy. In section III we establish the scale-locality of transfer in this 
decomposition. Section III. A establishes rigorous upper bounds on non-local transfers. 
Section III.B discusses decorrelation effects that lead to reduced non-local contributions 
to mean flux. Finally, section IV states our conclusions and two appendices provide some 
mathematical details. 

II Band-Pass Filtering via Smooth Coarse-Graining 
A The Low-Pass Filter Technique 

The filtering approach by smooth coarse-graining of the velocity field is standard in the 
LES literature. See, for example, papers'^^ and references therein. Any real-valued 
function G{r) may be chosen as a filter kernel, as long as it is sufficiently smooth, decays 
sufficiently rapidly for large r, and is normalized so that J d'^r G{v) = 1. It is assumed 
furthermore that G satisfies J d'^rrG{r) = and J d'^r \r\'^G{r) = 0{1), so that the 
main support of G is in the ball of radius 1 about the origin. Its dilation in d dimensions 
Ge{r) = i^'^G{r/i) shares these properties, except that it is supported mainly in the ball 
of radius i. One then defines the "coarse-grained" or (low-pass) filtered velocity field at 
length-scale i by 

u,(x) = Jdr G',(r)u(x + r). (1) 

In order to interpret this as a local space- average, G must also be positive, or G{r) > 
for all r. It is often convenient to choose G satisfying other desirable properties, such as 
radial symmetry. Note, for example, that G{r) may be chosen so that both it and its 
Fourier transform G(k) are positive and infinitely differentiable, with also G(k) compactly 
supported inside the ball of radius 1 in Fourier space and with G{r) decaying faster than 
any power as r oo. See Appendix 1. 

By applying the filtering operation to the incompressible Navier-Stokes equation, one 
obtains an evolution equation for the coarse-grained velocity u^: 

dtUi + (ui ■ V)ue = -Vvi - V ■ rf(u, u) + A u^, V ■ = (2) 

This is the same as the Navier-Stokes equation itself but with an additional contribution 
from the turbulent (Reynolds) stress, 

ri{u, u) = (uu)^ - UiUi. (3) 



4 



In the rest of the paper, we shall take the liberty of dropping the £ sub-script whenever 
there is no risk of ambiguity and use the shorthand r = r(u, u) to denote the stress ([3]). 
An energy budget for the large scales may be easily derived from ([2]), as: 

= -n-z/|vn|2 (4) 

where 

n = -idjUi)Tij (5) 

is what is usually called the sub-grid scale (SGS) dissipation or the SGS energy flux. It 
acts as a sink term in Q, representing the energy transferred from scales larger than i 
to the smaller (sub-grid) scales at point x in the flow. The terms inside the divergence 
represent energy transported in space. 

A small-scale energy or "subgrid kinetic energy" may be defined as 

^ = ^^(u; U) = ^T{Ui, Ui). (6) 

It is a positive quantity at every point in the flow if and only if the filtering kernel G{r) 
is positive for all r, as was proved by Vreman et al. (1994)^. Indeed, it can be rewritten 
as / dr G'(r)||u(x-|-r) — u(x)p, which is the energy density averaged over a region of size 
£ around x in a frame co-moving with the local large-scale velocity u(x)'I^. Furthermore, 
integrating ^f{ui,Ui) in space gives / (ix ||u(x)p — f ^|u(x)p, which is the total 
energy less the energy at large scales. It is easy to derive the energy budget of the small 
scales as: 

= U—uT^diUj, diUj). 

(7) 

For example, see Germano (1992)'^. The energy flux 11 now acts as source, representing 
the energy gained by the small scales from scales larger than £ at point x. 

B The Smooth Band-Pass Approach 

Following the same ideas, we may define a band-passed kinetic energy density at length- 
scales between £ and where £ > I, as 



^2' ' 



+ 9j {^\\l\^+p)Uj + UiT^j - Z/9j(^|u| 



dt-T{ui,Ui)+dj 



]rT{Ui, Ui)Uj+T{p, Uj) + \T{Ui, Ui, Uj)-udj]-T{Ui, Ui 



for filter functions G = Gj and G = Gj. The band-passed energy (8) is a straightforward 
generahzation of the small-scale kinetic energy presented in the previous sub-section. 
Positivity again holds, |r(u;u)(x) > 0, locally for all x when G'(r) > 0. Globally, this 
band-pass energy is 

j d'^x ^r(u; ,u) = j d'^x ^|u(x)|2 - J d^x ^|n(x)p 

^ ^ d'^k [1- |G(k£)p]|G(k£)|2j|u(k)|2, 



{2'kY j ' I V /I Ji V /I 2' 

by the Plancherel identity. This is clearly a reasonable representation of the energy in the 
wavenumber range [!/£,!/£]. If a sharp-spectral filter with G(k) = 9(2Tc — k) is used, then 
this band-pass energy reduces to I{2n/e<\k\<2TT/e}^'^^ which is the quantity 

considered in recent numerical studies^ ^ * ^ * ^ * ^ * ^ l 

One of the great advantages of the present approach is that very simple dynamical 
equations can be written for the band-pass energy densities above, which resolve the 
relevant physical process both in space and in scale. It is straightforward to show that 



+ ^T{Ui,Ui,Uk) +T(T{Ui,Uk),Ui) - udk^T(Ui,Ui 

= -SijT{ui, Uj) + T{ui, Uj)Sij - uT{diUj, diUj). (9) 

If the total energy (1/2) j |up remains finite in the limit ^ (as expected), then the 
viscous terms are easily seen to be negligible in a fixed band [£, for small v. The terms 
inside the space gradient once again represent space-transport of energy. The non- viscous 
terms on the righthand side of ([9]) represent the nonlinear transfer into and out of the 
band. Note that double overlining (■) is associated to the convolved filter G = G * G, 
which has a length-scale i ^ i foil £. Thus, the transfer terms have the meaning of a 
flux n into scales smaller than i and a reverse flux (—11) out of scales smaller than I. 

This is a good point to remark that the flux terms in eqs. (g, 0, and (g have 
all been deflned in a Galilean-invariant way, due to the subtracted mean terms in the 
definition ^ of the stress. Other definitions of an "energy flux" are possible, such as the 
"unsubtracted flux" 

^uns ^ —Q-Ui UiUj = n — (9j(-Kj|u|^). 

This differs from the standard SGS flux 11 by a total gradient, which can be included 
in the space-transport term. This "unsubtracted flux" is, in fact, often employed in 
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literature that considers the sharp-spectral filter. Unfortunately, this definition is not 
pointwise Galilean-invariant, so that the amount of "energy cascade" at any point in the 
fluid according to this definition would differ for observers moving at different uniform 
velocities! Kraichnan'^S', Spezialel^ and Germano'^ have all emphasized the importance 
of Galilean invariance. Our definitions are the unique ones which preserve the pointwise 
Galilean-invariance of energy flux. There are non-Galilean-invariant terms in our balance 
equations, but, as is physically natural, they are all associated to space-transport of 
kinetic energy. 

The previous methods may be used to introduce a band-pass decomposition of the 
energy at a geometric sequence of scales = p~"L with p > 1. The result is easiest to 
state for the special case of filter kernels G which satisfy the "S-type" , or sharp-spectral- 
like, condition that 

G(k) = li (lo) 



As we show in Appendix 1, it is possible to construct kernels of this type for which 
G(k) is C°° and thus G{r) decays faster than any power of r as r — > oo. We introduce 
Gn{r) = £~'^G'(r/£„) and the corresponding low-pass filter fn = Gn*f for any space 
function /(x). It is then possible to show that 

N 

\ / '"'^ " ^ / + T^«-i(un;u„) + ^ y r;v(u;u). (11) 

for any integer iV > 1. For the proof of this formula, and also a more general result 



without condition ( [10| ), see Appendix 2. It is easy to see that series (11) converges as 
N ^ oo whenever (1/2) J |up < oo. 

Using (|4]),([9]),([7]) one can infer evolution equations for the total band energies: 

(d/dt) j ^|uop = - /no, 

{d/dt) j ^x„_i(u„;u„) = j n„_i - j n„ for = 1, 

{d/dt) j ^rAr(u;u) = j li^ - p j TN{uij,Uij). 

Note that viscous terms may be neglected as z/ — in all of these equations except the 
final one, for n = N, where it balances the flux from the larger scales. If the turbulence is 
driven by a body force at length-scales > L, then an additional source term would appear 
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in the first equation for n = 0. Tliis liierarcliy makes apparent tlie stepwise nature of the 
cascade, if each of the flux terms n„ = — Vu„:rn(u, u) for inertial-range length-scales 
in depends only upon velocity modes of comparable scale. 

In a recent paper, Cheskidov et al. (2008)'^ have studied turbulent cascade dynamics 
by a band-pass decomposition of the velocity field with graded filters similar to ours. 
They derived thereby a number of estimates related to scale-locality of energy transfer. 
However, as noted by Domaradzki & Carati (2007)'^ the naive decomposition of the 
kinetic energy following from such a decomposition of the velocity field in general will 
contain off-diagonal terms: 



ur = / iuoi^ + ^ 5Z / 



1 /"l |2 1 

2/ ' ' 2 / ' 2 



n,m=l 

In the approach of Cheskidov et al. (2008)'^the summation is limited to indices with \n — 
m\ < 1 but still retains such off-diagonal terms of indefinite sign. It is thus not clear in this 
approach how precisely to identify the kinetic energy at a given length-scale. By contrast. 



our decomposition ( 11 ) of kinetic energy contains only non-negative, diagonal terms. The 
complications with the naive band-pass decomposition lead to related difficulties in the 
analysis of scale transfers of energy. For example, the "band-fiux" defined in eq.(27) of 
Cheskidov et al. (2008)^ is not the difference of the usual energy fluxes into and out of 
the band, as we found above in our approach, but contains an additional term [see their 
eq.(30)] with unclear physical interpretation, which those authors had to estimate. 



Ill Locality of Energy Cascade 

One of the advantages of our approach to band-pass decomposition of kinetic energy is 
that the proofs of scale-locality of Eyink (2005)'^^ directly apply. In fact, the flux quantities 
n„ were exactly the objects shown in that work to be local! To make connection with 
standard spectral approaches, we define the total transfer into the nth band by 

T„(u, u, u) = n„_i(u, u, u) - n„(u, u, u). 

Clearly, this quantity must be dominated by the contributions of local triads if this is true 
separately for both n„_i and n„. The total spectral transfer is often further decomposed 
as T„ = X]pTn,p = J2mp'^n,m,p into triadic contributions from distinct wavenumber 
bands. If we do so in the present approach of smooth coarse-graining, then it follows 
that transfer for non-local modal triads is negligibly small. 
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It should be emphasized that there are non-local triadic interactions in turbulent 
dynamics. For example, non-local interactions dominate in the space-transport of energy 
by convective sweeping. However, space-transport effects disappear upon integration 
over space and make no net contribution to the transfer of energy between scales. On 
the other hand, the energy fluxes n„ which are responsible for the inter-scale exchange 
of energy are scale-local, under well-defined conditions that exist in the inertial-range 
of turbulent flows. We give below a brief summary of the demonstration from Eyink 
(2005)111. 

A The Proof of Scale-Locality 

The energy flux 11^ is a triplet quantity that depends upon three velocity modes. This 
dependence may be made explicit by writing 

n£(u,v,w) = -Vu£ : r£(v,w). 

(In the discussion of locality it is important to indicate the scale £, which we do in this 
whole section). We say that the energy flux II^ = nf(u, u, u) is infrared local if the 
contribution from u at scales much larger than ^ is negligible. In other words, if u in 
the three arguments (u, v, w), or any subset thereof, were to be replaced by ua, then the 
result would be -C 11^ for /S.^ I. We can define a corresponding concept by employing 
the "small-scale" or high-pass filtered field 

U5 = u - U5, (12) 

which contains no modes at scales > 5. We say that II^ = n£(u, u, u) is ultraviolet local 
if replacing u by u'^ in the three arguments (u, v, w), or any subset thereof, were to give 
a result ^ II^ whenever 5 ^ This means that the contribution to energy flux from 
modes at scales much smaller than I is negligible. If a quantity is both infrared and 
ultraviolet local, then we say that it is (scale-)local. 

They key to proving scale-locality of the energy flux is the observation that all of the 
quantities Vu£,rf(u, u), and u'^ can be written entirely in terms of velocity-increments: 

5u(x; r) = u(x + r) — u(x) (13) 

for separation distances |r| < ^ (or some moderate multiple of I) and do not depend upon 
the absolute velocity u(x) itself. Heuristically, 

Vu^ ~ 0{5u{l)/l) (14) 
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t^(u,u)~0(5m2(£)) (15) 

~ 0{5u{l)) (16) 

The symbol ~ O stands for "same order-of-magnitude as." For rigorous details, see 
Eyink (2005)'i]. 

It thus becomes sufficient to show that velocity-increments themselves are scale-local. 
This is not generally true, for arbitrary solutions of the Navier-Stokes equations, but it is 
true under conditions that exist in the inertial-range of turbulent flows. The velocity fields 
in such turbulent solutions are not space-differentiable but only Holder continuous^. A 
field u is said to be Holder continuous at x with exponent < /i < 1 if its increments 
satisfy a rigorous big-0 upper bound 

|5u(x; r)| = O(r^). (17) 

The (maximal) Holder exponent at point x is the largest value h such that the inequality 



(17) still holds. We now observe that infrared contribution to the increment, (5ua(x; r) = 
ua(x + r) — ua(x), is negligible ii h < 1. Since the low-pass filtered field ua is very 
smooth, we can approximate it well by the leading term in its Taylor expansion: 

5nA(r; x) ^ r-VuA(x) ~ O = O(rA'^-i). (18) 

Compared with the full increment 5u(r) ~ 0{r^), the large-scale contribution is smaller 
by a factor of 0((r/A)^^^), which is indeed negligible if A ^ r and h < 1. We next 
observe that the ultraviolet contribution to the increment 5u'5(x; r) = u^(x + r) — u^(x) 



is negligible if h > 0. This is even easier to see, since, from (16), the high-pass filtered 
field itself is small: u'^ ~ 0{6u{6)). In that case, 

r) ~ 0{5u{5)) = 0(6''). (19) 

Compared with the full increment, the small-scale contribution is smaller by a factor of 
0{{S/r)^), which is also negligible if 5 <^ r and h > 0. All of these results hold pointwise, 
locally in space'^. 



One difficulty with the above pointwise analysis is that inequality (17) does not hold 



scaling lawElESI. The upper bound in \\18\j is vanishing for A ^ oo and the bound in 



(19) is vanishing for 6^0. However, because Holder continuity does not correspond to 
pointwise scaling, one cannot make precise statements about the fractional contributions 
from scales > A or < 5. On the other hand, it is well-known that there is scaling of 
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pth structure functions of velocity increments in the turbulent inertial-range. That is, 
(|5u(r)|P) ~ (9(M^^^(r/L)^p) where (■) = (l/V) jy{-) denotes volume- average. A global 



analogue of (17) thus holds for the Lp-norm || • ||p = ( | • 1^)^^^ as a true scaling law, i.e 



for some dimensionless constant Ap 



\\6n{r)\\pr^UrmsApir/Lr'' (20) 

with cTp = Cp/p. Experiments and simulations indicate that < dp < 1, and, in particular, 
(T3 = 1/3. The relation (20) corresponds to Besov regularitip^^ in place of the Holder 



continuity relation (17). It is not difficult to derive global Lp-mean analogues of IR- 



locality of velocity-increments 

||<5uA(r)||p = 0(rA'^^-^). (21) 
and UV-locality of velocity-increments 

||K(r)||p = 0(r^). (22) 



For details, see Eyink (2005)^. Together with (20) we may make the precise statement 
that the modes > A give a fractional contribution 0((£/A)^^°'p) and the modes < 6 a. 
fractional contribution 0{{6/iYp). 

These results for velocity-increments imply the locality of the energy flux. Here we 
shall just state two concrete estimates which correspond closely to classical statements 
of locality in the literature. First we consider the contribution to the flux from a large- 
scale strain, or Hf(uA,u, u) = — Vua/ : Te{u,u). Using the Holder inequality \{fgh)\ < 
ll/llalls'llsll^lla it ^^^J show from the foregoing that 

(H,(uA,u,u)) = 0(£(£/A)i-'^^). (23) 

We have used the relation e ~ uf^g/L in order to introduce the mean flux e. Next 
we consider the contribution from the small-scale stress, or H£(u, u'^, u'^) = — Vu^ : 
Te{u'g, u'^). Note that both of the modes in the stress must be small-scale together, if the 
fllter kernel G is "S-type" , because of wavenumber conservation constraints. Using again 
the Holder inequality and the previous estimates on increments, one finds that 

mn,n's,n's)) = 0iei5/ir^) (24) 



In deriving (23), (24) it was assumed that = 1/3, so that the exponents in both of the 



bounds are very close to 2/3. Also, to conclude that the non-local fiux contributions in 
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(23) and (24) are much smaller than (H^) for A ^ ^ and (5 <^ £, it must be assumed 
that the mean flux is non-zero. Of course, this is a completely realistic assumption in a 
constant-flux inertial range. 



B Decorrelation Effects 



At first sight, the results stated in (23) and (24) contradict the predictions of Kraich- 
nan'^El using closure calculations, which would yield a larger exponent 4/3 instead of 
2/3. However, one must keep in mind that our estimates above are only upper bounds 
or so-called big-0 bounds, which crudely replace all the velocity-increments with their 
absolute magnitudes. On the other hand, Kraichnan's estimates are more refined (but 
non-rigorous) scaling predictions which take into account cancellations in the signed con- 
tributions. See the heuristic discussion in Section 1 of Kraichnan (1971)'^ and the similar 
discussion by Tennekes and Lumley'^, sections 3.2 and 8.2. Their considerations can be 
carried over to our framework and then yield similar results. We now discuss this briefly. 
The easiest case to consider is infrared locality of the flux. In the average over space 

(n£(UA,U, u)) = I d'^xVuA/ ■ Ti{u,u), 

^ Jv 

the velocity-gradient Vua/ varies on a long length-scale ~ A whereas the stress T£(u, u) 
varies on the short scale ~ I. The average is only nonzero because of correlations between 
these two quantities, when the turbulence is either homogeneous or isotropic so that 
(Vu) = 0. An estimate of the correlation factor proposed inlSElis that 

p(VuA,r,) = 0(5(A)/5(£)) 

where S{t) ~ 0{5u{t) / 1) is the strain-rate at scale £. This can be plausibly justified by 
appealing to ergodicity to rewrite the space-averages as time-averages and then arguing 
that the correlation should be proportional to the ratio of the turnover rates at the two 
length-scales. Taking decorrelation into account yields an improved estimat^ 



The latter expression is a "multiscale structure-function" of a type much discussed in 
the literature. It can be estimated by means of a multiplicative cascade ansatz'^SES or, 

'^-'We remark here that the additional decorrelation resulting from space averaging of signed quantities 
can explain why Benzi et al.'^see a faster than expected decay rate of the multi-scale structure functions 
shown in their Figure 8. 
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equivalently, by so-called "fusion rules" '^^SSEU xhe result is that 



This gives an estimate for the flux contribution from the large-scale strain 



(n,(uA,U,u))~C> £ 



A 



if we take Cs = 1- The above estimate incorporates possible intermittency effects, but, 
assuming K41 scaling so that Ci = 1/3, we recover Kraichnan's predicted 4/3 exponent'^El. 

A similar discussion may be given for ultraviolet locality of the flux. Note first that 
Tiiyi!^, u'^) = (u'^u'^)^ for small enough 5, if the filter kernel G is "S-type" so that (u'^)^ = 
for 6 < i/ p. Thus, 

{n,(u, u^, u^)) = j^d'^x Vn,,, : u'^u'^ 

Since the velocity-gradient Vu^^^ varies on a long length-scale ~ i whereas the factors 
u'^ vary on the short scale ~ 6, we can expect decorrelation by a factor S{i)/S{6). This 
gives 

(n,(u,u'„u',)) ~ (^5«2(5)p(Vu,,T,)) ~ ^^{6u\i)6ui6)). 

Using the same heuristic estimation as before by the cascade ansatz or fusion-rule formula, 
one obtains the similar result 



(n,(u,u'„u',))~oU 



l+C: 



This again yields the 4/3 scaling predicted by Kraichnan^^S if we assume the K41 value 
Ci = 1/3. 

In fact, we can expect strong pointwise decorrelation effects in the ultraviolet contri- 
butions to the energy flux, due to local averaging over the scale i even without the global 
average over space. If the statistics of u'^ were isotropic in the ball of radius ~ i around 
point X, then (u'^u'^)^ = l(|u^P)^I. In reality, the statistics will not be exactly isotropic, 
because there is a local mean velocity-gradient Vu£(x) which distorts the smaller scales. 
Thus, we expect that 

, . 1- 
d' 



.s^sJe--AW)^^ + 0{S{i)/S{6)\u's\') 
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based on the same plausible reasoning as Tennekes & Lumlej^ and Kraichnan (1971)'^. 



The leading term does not make a contribution to energy flux, because of the zero trace 
of the velocity-gradient. Thus, pointwise, we expect that 



where n£(u, u, u) = O {5u^ {tj / C^) and we estimate 5u{5) = 0{6^) with local Holder ex- 
ponent h at X. This is considerably smaller than the rigorous pointwise upper bound 
estabhshed by Eyink (2005)'^. 

C Numerical Results 

The arguments in the preceding section are plausible but non-rigorous. It is therefore 
worthwhile to investigate the predicted decorrelation effects in numerical simulations. 
We present here the results of a 512^ DNS of Navier-Stokes turbulence, whose details 
are discussed in paper II. We use Gaussian filters Ge{r) = (vr/f^)^/^ exp(— vr^r^/^^) to 
scale-decompose the simulated velocity field. To facilitate comparison with results for 
spectral filters we employ an "equivalent wavenumber" K with = 2ti/K. 

We first check the decorrelation argument in the case of IR locality, by measuring both 
the mean absolute flux (|n£^(u£Q, u, u)|) and the mean signed flux (n£^(u£Q, u, u)), for 

> 2£j^. The former is immune to cancellation effects, whereas the second should 
experience decorrelation from space averaging. The plots in figure [T] are consistent with 
a Q^/^ scaling for absolute flux and a Q^^^ scaling for signed flux, exhibiting the faster 
decay due to decorrelation effects. The fitted IR exponents are slightly smaller than those 
predicted, consistent with the results of Domaradzki & Carati^^. 

We also check the decorrelation argument for UV locality by measuring the abso- 
lute flux (|Vu^^ £^- : u^^u^^l) and the signed flux (Vu^^^^^ : u^^u^^), for Ip < ixf^- 
The plots in figure [2] are again consistent with respective decays of P~^/^ and P~^/^, 
as expected from our arguments. The fitted UV exponents are somewhat larger than 
predicted, also as observed by Domaradzki & Carati'^. 
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Figure 1: For K = 100, we plot as functions of Q < K/2 both (|n£^(u£^, u, u)|)/ (in^^^ |) 
(o) and {IIi^{uiq,u,u)) / {Uij^) (x) using a gaussian filter. Straight lines of slopes 2/3 
and 4/3 are included for reference. They extend over the fitting region which gives a slope 
0.7 for the absolute flux and 1.1 for the signed flux. The difference in slope demonstrates 
the decorrelation effect due to space averaging. 
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Figure 2: For K = A, we plot as functions of P > 2K both (jVu^^^^^ : u^^u^^|)/(|n£^|) 
(o) and (Vu^^^£^ : u'^^u'^^) / {Uij^) (x) using a gaussian filter. Straight lines with slopes 
2/3 and 4/3 are added for reference. They extend over the fitting region which gives a 
slope 0.9 for the absolute flux and 1.9 for the signed flux, demonstrating the decorrelation 
effect due to space averaging. The kink at very large P is due to negative values of the 
function when plotted on a log-log graph. 
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IV Conclusions 



Let us summarize our main results. First, we have shown how to obtain a band-pass de- 
composition of kinetic energy with a smooth, graded filter, corresponding to length-scales 
in a geometric sequence. This decomposition permits simultaneously resolution of the 
turbulence dynamics both in scale and in physical space. Second, we have demonstrated 
that the interband transfers in this decomposition are dominated by local triads of modes, 
assuming only the scaling properties of velocity-increments that are expected (and ob- 
served) to hold in turbulent flows. Finally, we have explained the apparent discrepancy 
between our rigorous upper bounds and the scaling predictions of Kraichnan'^El, as due to 
cancellations that appear in averages over space. These decorrelation effects have been 
confirmed by results from a DNS of homogeneous, isotropic Navier-Stokes turbulence. 

We should emphasize the importance of the pointwise results on scale locality that 
were deduced here and in the previous workf^. If local triadic interactions dominated the 
energy cascade only "on average" and not pointwise in space, for individual flow realiza- 
tions, then there would be scant basis for the usual hypothesis of small-scale universality. 
Higher-order statistics in particular could show non-universal effects from such non-local 
interactions, if they were sizable. Fortunately, this is not the case. It should be noted, 
however, that simultaneous resolution of the turbulence fields both in space and in scale 
requires that our decomposition employ bands of constant width on a logarithmic scale. 
Our upper bounds on non-local contributions are exponentially small in the logarithmic 
band-number. However, this corresponds to only a slow, power-law decay in length-scale 
ratios, in agreement with earlier conclusions of Kraichnanl^ZE Tennekes & Lumley'^ 
about the "diffuseness" or "leakiness" of the energy cascade. Thus, very high Reynolds 
numbers and long inertial ranges could be required before universal statistics are obtained 
at small scales. 

A final conclusion of this paper is that smooth band-decompositions which provide 
space-scale resolution have energy transfer between bands due mainly to local triadic 
interactions. It remains to account for the disagreement with DNS studies^ ° * ^^ * ^^ * ^^ * ^ ^ 
that claim to observe transfer predominantly by non-local triads. This is the subject of 
the following Part II of this work. 
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for assistance in developing the simulation code used in this work. Computer time was provided 
by the Digital Laboratory for Multi-Scale Science at the Johns Hopkins University. This work 
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Appendix 1: Filter Kernels with Compact Spectral 

Support 



In this appendix we explain briefly some elementary methods to construct filter kernels 
with compact spectral support, used to simplify several arguments in the text. 

The first construction yields a filter kernel whose Fourier transform G is C°° and 
supported inside the ball in Fourier space with |k| < 1. Thus, G is decaying faster than 
any power of |r| as |r| oo and is analytic in r. We shall show that both G and G 
may be chosen, furthermore, to be nonnegative and rotationally symmetric. To begin 
the construction, define a function G'o(k) of k = |k| which is C°° and supported inside 
the ball |k| < 1/2. This is a standard problem in analysis. For example, one may take 



I |k| > 1/2. 

To show that this function is infinitely differentiable it is enough to check that radial 
/c-derivatives of all orders are zero at k = 1/2, both from the left and the right. The 
above property implies that the inverse Fourier transform 



Go(r) = ^ j rf"A;Go(k)cos(k.r) 



decays faster than any power \r\~^ as |r| oo. Furthermore, Go{r) is real- analytic, 
i.e. has a convergent power series representation with an infinite radius of convergence. 
The function Go{r) is not everywhere positive, but G{r) = [Go{r)]'^ is so. Also, G{r) is 
real-analytic and decays faster than any power as |r| ^ oo. It follows easily that 



1 ' 



G(k) = — J dVGo(p)Go(k-p) 

is positive, C°° and compactly supported in the ball |k| < 1. In Figsjs] and |4] we have 
plotted the 1-dimensional versions of G and G that follow from this construction, with 
the choice of Gq given above. 

As discussed in the text, it is often convenient to use filter kernels of "S-type", which 



satisfy the condition (10) in the text for some p > 1. To construct such a kernel, we first 
define a function H[k) in Fourier space by 

< A; < 1 



Hik) = { exp(^-(^3^^j l<k<p 
k>p 
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Figure 3: A filter kernel G{x) which is positive, real analytic, and decays faster than any 
power at large distances. 



G 




Figure 4: The Fourier transform G{k) of the kernel in Figjsjis positive, C°°, and com- 
pactly supported in the unit ball. 
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which is C°° and compactly supported in the interval [1, p]. We can now define G{k) by 

G{k)^l-^j\pHip) 
with n = J[dpH{p) > 0. Clearly, 

^ ^ \ iik> p 

G(k) = G{k) is also C°° so that its d-dimensional inverse Fourier transform G{r) decays 
faster than any power as r ^ oo. Furthermore, G{r) is real-analytic in r, because 
G{k) is compactly supported. 

Note, however, that filter kernels G of "S-type" cannot be non-negative, since 



J d'^r\r\^''G{r) = (^(-A)"g) (k) 



= 

k=0 



for all positive integers n. In that respect, they are similar to the sharp spectral filter, 
whose physical-space kernel is proportional to a Bessel function that also takes on negative 
values. The primary advantage of the above construction is that the kernel G is much 
better localized in physical space than is the kernel of the sharp-spectral filter, which 
decays only as a power-law for r — > oo. 
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Appendix 2: Band-Pass Decomposition of Quadratic 

Integrals 



We show here how to construct a band-decomposition of quadratic integrals for a graded 
filter. Let us first make some definitions. We denote the filter function at length-scale 

in = p-'^i by 

for a smooth filter kernel G in space dimension d. The corresponding low-pass filter is 



/„(x) = (G„*/)(x) = J t/rG„(r)/(x + r) 



for any space function /(x). In our development below we need also the repeated low-pass 
filter 

7m,n(x) = {Gn * Gn+1 *■■■* Gm-l * Gm * /)(x) 

for any n < m. Note that we have successively removed larger and larger length-scales, 
so that this quantity is another type of low-pass filter at length-scale for any n < m. 
Finally, we define 

^n(/,^) = {fg)n-ln9n- 

This is the Germano (1992) "generalized central moment" of 2nd-order, at length-scale 

'"HI- 

We observe the following identity for arbitrary space functions /, g and integer iV > 0: 

N 

if9)N,o ~ f n,q9n,o + 

n=l 

It is not hard to see that successive terms cancel in this telescoping sum. This iden- 
tity generalizes the well-known Germano identity'^, to which it reduces when = 1 : 
'Ti,o{f,g) = To{fiy9i) + {'Ti{fy9))o- Integration over space gives an exact band decompo- 
sition of a general quadratic integral, 

N 

f9= 7n,09n,0 + Y1 / '^n~^(fN,n,9N,n)+ / ^Jv(/,^)- (A-1) 

Here we have used the fact that / / = / / for any filtered quantity. The physical 



interpretation of (A-1) is that the first term corresponds to the contribution from length- 



scales > io, the next terms for n = 1,...,N correspond, respectively, to contributions 
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from length-scales between in-i and in, and the final term corresponds to contributions 
from length-scales < i^- 

A simplified formula holds for filter kernels of "S-type". In this case, /^,^ = /„ for 
all n < N. Therefore, 

N 

J J J J 



This gives the formula (11) of the text if one takes f = g = Ui and then sums over i. 
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